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Abstract 

A  study  has  been  conducted  on  the  effects  of  flow  distributor  and  diffusion  layer  geometries  and  morphologies  on  the  transport  phenomena 
of  reactant  gas  from  flow  channels  to  gas  diffuser,  and  on  the  performance  of  a  proton  exchange  membrane  fuel  cell.  The  emphasis  of  the 
analysis  is  placed  on  the  effects  of  the  cross-section  of  the  channel  of  the  anode  flow  distributor  and  the  porosity  of  the  gas-diffusion  layer. 
By  applying  a  two-dimensional  mass  transport  model,  the  effects  of  the  channel  width  fraction,  X  =  ic/ib,  the  channel  number,  N ,  the 
porosity  of  the  diffuser  layer,  s,  and  the  surface  overpotential  of  the  catalyst  layer,  p,  on  the  resultant  current  density  are  investigated.  The  cell 
performance  is  evaluated  by  the  predicted  voltage-current  density  curves.  It  is  disclosed  that  an  increase  in  either  X,  Nor  s  may  lead  to  a  better 
cell  performance.  At  relatively  low  overpotential,  better  uniformity  in  current  density  distribution  along  the  width  of  the  cell  (v-direction) 
can  be  attained.  Based  on  the  present  results,  a  correlation  of  overpotentials  is  proposed,  namely,  pCY  =  o.1652.0  0935A°-2033£:0  0856,  which 
may  serve  as  a  guideline  to  justify  the  validity  of  a  simplified  one-dimensional  model. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  is 
a  most  promising  energy  source  due  to  its  environmen¬ 
tal  attributes  as  well  as  its  high  efficiency  The  latter  fea¬ 
ture  depends  strongly  on  the  operating  conditions,  transport 
phenomena  in  the  cells,  electrochemical  reaction  kinetics, 
mechanical  design  and  manufacturing  process,  etc.  Crucial 
among  these  is  the  reactant  gas  transport  in  the  flow  dis¬ 
tributor  and  the  gas-diffusion  layer  before  the  reactant  gas 
reaches  the  catalyst  layer  to  carry  out  reaction.  Different 
types  of  flow  channel,  such  as  the  serpentine  channel,  multi¬ 
ple  channels  in  parallel  type  and  interdigited  channels,  have 
been  analyzed  in  order  to  reach  optimization  of  the  cell  per¬ 
formance.  The  geometric  parameters  and  the  morphologies 
of  the  gas-diffusion  layer  and  the  catalyst  layer  have  also 
been  investigated. 


*  Corresponding  author.  Fax:  +886-4-24516246. 

E-mail  address:  cysoong@fcu.edu.tw  (C.Y.  Soong). 

0378-7753/$  -  see  front  matter  ©  2003  Elsevier  B.V.  All  rights  reserved. 
doi:10.1016/j.jpowsour.2003.07.017 


West  and  Fuller  [1]  studied  the  effects  of  the  rib  spac¬ 
ing  of  bipolar  plates  and  the  thickness  of  the  gas-diffusion 
electrode  on  cell  performance.  It  was  found  that  the 
two-dimensional  effect  only  slightly  alters  the  half-cell 
potential  for  a  given  applied  current,  but  has  a  significant 
influence  on  water  management.  Jordan  et  al.  [2]  examined 
the  influence  of  diffusion-layer  morphology  on  cell  perfor¬ 
mance.  A  model  of  the  hydrophobilicity  and  the  porosity 
of  the  gas-diffusion  layer  was  used  to  explain  the  influence 
of  the  diffusion-layer  morphology.  In  other  work  [3],  these 
authors  obtained  results  which  showed  that  a  low-porosity 
acetylene  black  enhances  water  removal  from  the  membrane 
electrode  assembly  (ME A),  and  thus  gives  an  improvement 
in  gas-diffusion.  Guran  et  al.  [4]  proposed  a  mathematical 
model  to  obtain  an  analytic  solution  of  the  mass  transport 
of  reactant  gas  in  a  half-cell,  in  which  the  effects  of  the 
porosity  and  the  tortuosity  of  the  gas-diffusion  layer  and 
the  catalyst  layer  were  explored. 

To  obtain  details  of  the  flow  distribution  and  gas-diffusion, 
numerical  computations  have  been  carried  out  recently  by 
a  few  investigators.  For  example,  Dutta  et  al.  [5]  solved 
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Nomenclature 

D  Mass  diffusivity  (m2  s-1) 

F  Faraday  constant,  96,487  C  mol-1 

I  current  density  (Am-2) 

/o  maximum  current  density  at  v  =  0  (Am-2) 
j  transfer  current  density  (Am-2) 

t b  base  width  of  the  flow  distributor  (m) 
ic  gas  channel  width  (m) 

is  shoulder  width  (m) 

Mo2  molecular  mass  of  oxygen,  32.0  kg  mol-1 
n  number  of  electrons  in  transfer  current 

density  relation 

N  number  of  flow  channels 

P  pressure  (Pa) 

R  universal  gas  constant,  8.314  Jmol-1  K-1 

tc  thickness  of  catalyst  layer  (m) 

td  thickness  of  diffusion  layer  (m) 

tm  thickness  of  membrane  layer  (m) 

T  temperature  (K) 

V  cell  voltage 

w  oxygen  mass  fraction 

wq  oxygen  mass  fraction  at  inlet 

Greek  symbols 

a  transfer  coefficient  for  the  reaction 

£  \  porosity  of  diffusion  layer 

£2  porosity  of  catalyst  layer 

£i  volume  fraction  of  ionomer  in  the  catalyst  layer 
t)  surface  overpotential  (V) 

k  reaction  rate  constant  (mol  m-3  s-1) 

A  channel  width  fraction,  ijib 

Aw  membrane  water  content,  kmol  F^O/kmol  S03~ 
p  density  of  oxygen  (kg  m- 3 ) 

o  ionic  conductivity  of  the  ionomer  (Q~l  m-1) 
x  i  tortuosity  of  diffusion  layer 

T2  tortuosity  of  catalyst  layer 

Ti  tortuosity  of  the  ionomer  in  catalyst  layer 

<fi  ionomer  phase  potential  (V) 


three-dimensional  Navier-Stokes  equations  with  a  electro¬ 
chemical  relation  model  for  the  mass  source/sink  terms  to 
examine  the  mass  exchange  in  serpentine  flow  channels  at 
both  the  cathode  and  the  anode.  To  achieve  optimization  of 
the  flow  field  and  to  enhance  fuel  cell  performance,  Hon- 
tanon  et  al.  [6]  focused  attention  on  the  flow  distributor  rather 
than  the  whole  cell.  Two  types  of  flow  distributors,  were 
evaluated,  namely,  grooved  plate  and  porous  material,  by 
employing  Navier-Stokes  equations  as  the  flow  model  and 
Michaelis-Menten  type  two-step  kinetics  as  the  gas  reaction 
model  at  the  anode.  The  computational  results  showed  that 
fuel  consumption  can  be  increased  by  decreasing  the  per¬ 
meability  of  the  flow  distributor,  especially  when  the  per¬ 
meability  decreases  to  a  value  below  that  of  the  anode.  In 


terms  of  the  reactant  gas  utilization  based  on  the  same  pres¬ 
sure  drop,  it  was  claimed  that  porous  materials  are  more 
advantageous  than  grooved  plates.  For  flow  distribution  in 
the  channel  networks  of  typical  planar  fuel  cell  layers  and 
stacks,  Kee  et  al.  [7]  have  recently  proposed  a  generalized 
computational  model,  in  which  only  mass  and  momentum 
transfer  are  considered. 

The  objective  of  the  present  work  is  to  study  the  effects  of 
the  flow  distributor  and  the  diffusion  layer  morphologies  on 
the  transport  of  reactant  gas  from  the  flow  channels  to  the  gas 
diffuser,  and  on  the  performance  of  a  PEMFC.  Emphasis  is 
placed  on  the  effects  of  the  cross-section  of  the  flow  in  terms 
of  the  ratio  of  the  channel  width  i c  to  the  base  plate  width 
£b,  he.  A  =  ic/i b,  and  the  morphology  of  the  gas-diffusion 
layer  in  terms  of  the  porosity  of  the  layer.  By  applying  a 
mass  transport  model,  the  effects  of  these  two  parameters 
and  the  overpotential  on  the  resultant  current  density  are 
investigated  and  the  cell  performance  is  evaluated  by  the 
voltage-current  density  or  performance  curves.  Based  on  the 
results  with  various  overpotentials,  a  correlation  is  proposed 
to  serve  as  a  guideline  for  justification  of  the  validity  of  a 
simplified  one-dimensional  model. 


2.  Theory 

2.7.  Model  development 

A  schematic  diagram  of  the  cross-sectional  view  of  a 
half-PEMFC  is  shown  in  Fig.  1.  The  shoulder  area,  which 
is  essential  for  electronic  transport  from  the  electrode  to  the 
current-collector,  affects  the  fuel  cell  performance  by  chang¬ 
ing  the  ohmic  and  contact  resistance  in  the  cell.  Therefore, 
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Fig.  1.  Half-cell  model  of  a  PEMFC.  FC:  flow  channel;  CS:  channel  shoul¬ 
der;  CC:  current  collector;  BP:  bipolar  plate;  GDL:  gas-diffusion  layer; 
CL:  catalyst  layer;  M:  membrane;  B1-B7:  boundary  numbers. 
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it  is  important  to  investigate  the  effects  of  shoulder  area  on 
the  transport  phenomena  of  the  reactant  through  the  cathode 
and  on  the  ohmic  resistance  of  the  flow  distributor.  More¬ 
over,  the  area  ratio  between  flow  channel  and  shoulder  needs 
to  be  optimized  for  both  gas  transports  and  electronic  re¬ 
quirements. 

Oxygen  gas  transport  through  a  half-PEMFC  is  simu¬ 
lated.  Gas  transport  from  the  flow  channels  to  the  catalyst 
is  by  diffusion,  which  is  driven  by  the  concentration  gradi¬ 
ent.  A  steady-state,  two-dimensional  transport  equation  for 
the  oxygen  mass  fraction,  w ,  in  the  diffusion  layer  can  be 
expressed  by: 


9 

dx 


where  p  is  the  gas  density;  D  the  mass  diffusivity;  £\  and 
x  i  are  the  porosity  and  tortuosity  factor  of  the  diffusion 
catalyst  layer,  respectively.  In  the  catalyst  layer,  the  oxygen 
is  consumed  and  the  transport  equation  becomes: 


9 

dx 


27fM o2=°  (2) 


where  j  is  the  transfer  current  density;  n  the  number  of  elec¬ 
trons  in  the  transfer  current  density  relation;  F  the  Faraday 
constant;  £2  and  T2  are  the  porosity  and  the  tortuosity  factor 
of  the  catalyst  layer,  respectively. 

Assuming  that  the  oxygen  reduction  reaction  (ORR)  is 
irreversible  and  first  order  in  oxygen,  and  that  electroneu¬ 
trality  holds  in  any  representative  elementary  volume,  the 
cathode  transfer  current  density  can  be  expressed  as: 


j  —  uFkw  exp 


Substituting  Eq.  (4)  into  (2),  gives: 


2.2.  Boundary  conditions 


The  conditions  subjected  to  the  boundaries  numbered 
B1-B7  in  Fig.  1  are  described  as  follows.  The  boundary 
conditions  for  the  oxygen  mass  fraction,  w,  are: 


dw 


dy 

dw 

dy 


:  W0, 

at  boundary  B 1 

(8a) 

=  0, 

at  boundaries  B2  and  B5 

(8b) 

=  0, 

at  boundaries  B3  and  B7 

(8c) 

dw  dw 

W-  ~  u;+,  s\l  D —  =  sl2D  — ,  at  boundary  B4. 

dy  dy 


(8d) 


The  boundary  conditions  for  potential  0  are: 


90 

dy 

90 
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=  0,  at  boundary  B4 


=  0,  at  boundaries  B3  and  B7 
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3  dy  _  dy 

0  =  0,  at  boundary  B 6. 

2.3.  Evaluation  of  cell  potential 


,  at  boundary  B5  (9c) 


+ 


(9d) 


With  the  potentials  given  by  the  above  model,  the  cell 
potential  V  can  be  calculated  by  the  following  equation: 


v  =  Voc  -  r)  -  A0 


(10) 


The  expression  for  the  membrane  phase  potential  is: 


2  (a* )  +  2 
dx  \  dx  J  dy 


In  the  catalyst  layer  with  the  presence  of  transfer  current 
density,  the  phase  potential  equation  is: 


2  (e?a*t)  +  2  (e«a*) 

dx  \  dx  J  dy  \  dy  J 


where  £\  and  t\  are  the  volume  fraction  and  the  tortuosity 
factor  of  the  ionomer  in  the  catalyst  layer;  Aw  the  membrane 
water  content;  o  is  the  membrane  ionic  conductivity  and  is 
evaluated  by: 


a  =  (0.5139AW  —  0.326)  exp 


1268 


X  10“2 


where  Voc  is  the  ideal  cell  potential;  77  the  surface  overpo¬ 
tential;  A0  is  the  average  phase  potential  difference  and  can 
be  expressed  as: 


(11) 


in  which  lb  =  lc  +  ls  is  the  base  width  (one-half  of  a 
channel- shoulder  pitch),  and  04  and  06  denote  phase  poten¬ 
tials  at  boundaries  B4  and  B6,  respectively. 


3.  Results  and  discussion 

The  model  described  above  is  solved  by  a  finite  differ¬ 
ence  method  with  a  uniformly  distributed  grid  having  201 
and  221  points  in  the  v  and  y  directions,  respectively.  The 
values  of  the  parameters  considered  in  the  present  study  are 
listed  in  Table  1 .  The  numerical  results  are  presented  to  dis¬ 
cuss  the  effects  of  some  physical  and  morphological  param¬ 
eters.  The  geometry  is  characterized  by  a  parameter  named 
the  channel  width  fraction,  A  =  lj lb,  of  which  three  val¬ 
ues,  A  =  0.25,  0.5,  and  0.75  are  considered  to  explore  the 
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Table  1 


Parameters  used  in  present  analysis 


Parameters 

Notation 

Value 

Temperature  (K) 

T 

333 

Pressure  (Pa) 

P 

303975 

Reaction  rate  constant  (molm-3  s-1) 

K 

1.5  x  10“5 

Oxygen  mass  fraction  at  inlet 

wo 

0.17 

Gas  channel  width  (m) 

4 

(a)  4-channel,  5  x  10~4 

(b)  2-channel,  1.0  x  10-3 

(c)  Single-channel,  2.0  x  10-3 

Channel  shoulder  width  (m) 

4 

(d)  4-channel,  5  x  10-4 

(e)  2-channel,  1.0  x  10-3 

(f)  Single-channel,  2.0  x  10-3 

Thickness  of  diffusion  layer  (m) 

u 

'■t 

i 

o 

H 

X 

in 

Thickness  of  catalyst  layer  (m) 

tc 

5  x  10"5 

Thickness  of  membrane  (m) 

tm 

2.0  x  10“4 

Porosity  of  GDL 

£t 

0.8 

Tortuosity  of  GDL 

G 

1.5 

Porosity  of  catalyst  layer 

£2 

0.3 

Tortuosity  of  catalyst  layer 

G 

1.5 

Volume  fraction  of  ionomer  in  CL 

£3 

0.3 

Ionomer  tortuosity  in  CL 

G 

1.5 

Faraday  constant  (Cmol-1) 

F 

96487 

Universal  gas  constant  (kJkmol-1  K-1,  Jmol-1  K-1) 

R 

8.314 

Molecular  mass  of  oxygen  (kgkmol-1) 

M 

32.0 

Number  of  electrons  in  relation  participating 

n 

2 

Current  density  (Am-2) 

I 

Surface  overpotential  (V) 

7] 

Ionomer  phase  potential  (V) 

0 

geometric  effects  of  the  flow  channel  on  the  performance. 
In  addition,  the  number  of  flow  channels,  N ,  is  also  an  im¬ 
portant  parameter  worthy  of  investigation.  For  a  fixed  chan¬ 
nel  width  fraction,  A  =  0.5,  three  configurations  of  N  =  1, 
2,  and  4  are  studied.  Finally,  to  explore  the  effects  of  the 
porosity  of  the  gas-diffusion  layer,  various  values  of  poros¬ 
ity  are  considered.  Without  loss  of  generality,  the  porosities 
of  the  catalyst  and  membrane  have  the  same  value  as  that 
of  the  gas-diffusion  layer,  i.e.  £2  =  s\  =  £  1,  and  porosity  in 
the  range  of  0.3  <  £\  <  0.6  are  examined.  In  the  following 
discussion,  unless  otherwise  stated,  the  results  are  generated 
at  the  condition  s\  =0.3. 

3.1.  Effect  of  channel  width  fraction 

The  contours  of  the  oxygen  mass  fraction  in  the  diffusion 
layer  (0  <  y  <  5  x  10  4)  and  the  catalyst  layer  (5  x  10-4  < 
v  <  5.5  x  10-4)  of  s\  =  0.3  in  the  case  of  TV  =  2,  rj  =  0.2 
and  A  =  0.25,  0.5  and  0.75  are  presented  in  Fig.  2.  The  bold 
line  on  the  bottom  right  of  each  contour  marks  the  axial  lo¬ 
cation  of  the  shoulder  (or  rib)  of  the  channel.  In  the  region 
corresponding  to  the  channel  central  region  (left  portion  of 
the  contour),  a  relatively  higher  oxygen  mass  fraction  ap¬ 
pears  and,  then,  the  oxygen  diffuses  towards  the  porous  ma¬ 
terial  that  covers  the  rib  where  there  is  no  oxygen  source  on 
the  interface  y  =  0.  The  concentration  gradient  of  the  oxy¬ 
gen  is  more  noticeable  in  the  case  of  a  smaller  channel  width 
fraction,  i.e.  X  =  0.25,  and  with  increasing  channel  width, 


the  gradient  becomes  moderate.  The  phase  potential  distri¬ 
butions  in  the  catalyst  layer  (5  x  10  4  <y  <5. 5x10  4)and 
the  membrane  (5.5  x  1CT4  <  y  <  7.5  x  10  4)  of  porosity 
£  =  0.3  are  presented  in  Fig.  3.  It  is  found  that  an  increase 
in  the  channel  width  fraction  alleviates  the  potential  gradi¬ 
ent  in  the  v-direction.  Additionally,  at  a  large  value  of  A,  the 
total  difference  between  the  potentials  06  at  y  =  7.5  x  10-4 
and  04  at  y  =  5  x  10-4  along  the  v-direction,  A 0  = 
J 0  (06  —  04)  &x/l\y,  becomes  small  and,  from  Eq.  (10),  the 
cell  potential  increases.  This  phenomenon  is  attributed  to  a 
more  uniform  distribution  of  oxygen  (Fig.  3). 

The  current  density  distribution  at  the  interface  between 
bipolar  plate  and  gas-diffusion  layer,  i.e.  boundaries  B1  and 
B2  in  Fig.  1,  under  various  given  conditions  of  surface  over¬ 
potential,  rj ,  is  shown  in  Fig.  4.  With  periodic  conditions 
at  left  and  right  sides,  i.e.  boundaries  B3  and  B7,  the  ef¬ 
fects  of  relative  channel  width  are  explored  by  examining 
the  current  density  distributions  in  the  v-direction.  It  is  ob¬ 
served  that  the  distribution  of  the  current  density  in  the 
v-direction  at  the  interface  of  the  gas-diffusion  layer  and 
the  catalyst  layer  can  be  uniform  at  low  overpotentials.  At 
a  high  overpotential,  the  surface  reaction  rate  becomes  fast 
and  the  relatively  weak  diffusion  leads  to  a  non-uniform 
current  density.  This  situation  is  prevailing  over  the  various 
channel  width  fractions  considered,  i.e.  A  =  0.25,  0.5,  and 
0.75. 

Current  density  distribution  with  variation  of  overpoten¬ 
tial  is  plotted  in  Fig.  5.  The  data  clearly  show  the  level 
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Fig.  2.  Oxygen  mass  fraction  contours  in  diffusion  and  catalyst  layer  with  N  =  2,  rj  =  0.25,  and  (a)  7  =  0.25;  (b)  7  =  0.5;  (c)  7  =  0.75. 
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Fig.  3.  Ionomer  phase  potential  contours  in  catalyst  layer  and  membrane  with  N  =  2,  rj  =  0.25,  and  (a)  A  =  0.25;  (b)  A  =  0.5;  (c)  A  =  0.75. 
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Fig.  4.  Effects  of  surface  overpotential  on  current  density  distribution  along  x-direction  with  N  =  2  and  various  channel  width  fractions,  (a)  A  =  0.25; 
(b)  X  =  0.5;  (c)  X  =  0.75. 


and  the  non-uniformity  of  I  distribution  with  increasing 
overpotential.  The  maximum  current  density  which  appears 
at  the  location  v  =  0  is  denoted  by  /o.  The  ^-dependence 
of  /o  is  presented  in  Fig.  6.  It  is  found  that,  as  overpotential 
is  increased  to  enhance  the  surface  reaction,  the  maximum 
current  density  changes  dramatically.  Increase  in  the  over¬ 
potential  enhances  the  surface  reaction.  The  local  maximum 
current  density  /o  is  a  parameter  can  be  used  to  characterize 
the  level  of  the  current  density  in  the  fuel  cell  model.  The 
data  in  Fig.  6  disclose  that  the  variation  of  Iq  can  be  divided 
into  three  regimes: 

(i)  slow  reaction  regime  (rj  <  0.12) — in  which  the  current 
density  is  low  and  the  increase  in  77  has  little  effect  on 

/o; 

(ii)  high  sensitivity  regime  (0.12  <  rj  <  0.3) — the  regime 
is  characterized  by  remarkable  changes  in  Iq  with  in¬ 
creasing  K)\ 


(iii)  saturation  regime  (77  >  0.3) — where  Iq  approaches  an 
asymptotic  value  and  further  increase  in  77  has  little 
effect  on  Iq. 

The  correlations  of  current  density  with  overpotential  in 
Fig.  6  demonstrate  that  the  influences  of  the  channel  width 
fraction  on  Iq  are  just  moderate.  The  Iq  versus  77  curve 
at  A  =  0.5  almost  coincides  with  that  at  A  =  0.75.  This 
implies  that  the  effect  of  A  on  7o  vanishes  at  A  >  0.5.  Al¬ 
though  the  maximum  values  of  current  density  for  the  two 
geometries  are  extremely  close,  the  distributions  shown  in 
Fig.  4(b)  and  (c)  provide  evidence  of  better  performance 
with  A  =  0.75  than  with  A  =  0.5.  This  is  attributed  to 
the  higher  local  values  of  current  density  and  the  bet¬ 
ter  uniformity  attained  in  a  case  of  larger  channel  width 
fraction.  To  explore  the  influence  of  channel  width  on  per¬ 
formance,  the  I  versus  V  correlations  at  A  =  0.25,  0.5,  0.75 
and  1.00  are  presented  in  Fig.  7.  These  curves  demonstrate 
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Fig.  5.  Quantitative  representation  of  current  density  distributions  along  x-direction  corresponding  to  cases  in  Fig.  4.  (a)  A  =  0.25;  (b)  A  =  0.5;  (c)  A  =  0.75. 


that  the  larger  is  the  channel  width  fraction,  the  better  is 
the  performance.  Without  consideration  of  the  details  of  the 
electrochemical  reaction  and  other  factors,  the  results  are 
reasonable  from  a  geometric  viewpoint. 

3.2.  Effect  of  channel  number 

To  explore  the  influence  of  the  flow  channel  number  on 
the  diffusion  of  gas  reactants  and  the  performance  of  the  fuel 


cell,  the  models  with  a  fixed  channel  width  fraction,  X  =  0.5, 
are  used.  The  effect  of  overpotential  on  the  current  density 
distribution  with  channel  number  N  =  1 ,  2  and  4  is  shown  in 
Fig.  8,  from  which  the  effect  of  channel  number  can  also  be 
extracted.  From  a  comparison  of  the  local  distribution  of  the 
current  density  at  a  fixed  overpotential,  it  is  found  that  better 
uniformity  can  be  obtained  for  a  flow  channel  design  with 
a  larger  channel  number.  The  inherent  distributive  nature  of 
the  multi-channel  design  is  the  major  reason.  Three-regime 
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Fig.  6.  Correlation  of  local  maximum  current  density  and  surface  overpotential  with  various  channel  width  fractions. 


nature  of  the  maximum  current  density  curves,  i.e.  /o  versus 
0  plots,  defined  above  is  clearly  seen  in  Fig.  9.  The  effects  of 
channel  number  on  the  value  of  /o,  like  that  of  the  channel 
width  fraction,  have  little  significance. 

The  influence  of  the  channel  number  on  the  performance, 
i.e.  the  I  to  V  correlations  in  the  cases  of  N  =  1,2,  and  4 
are  presented  in  Fig.  10.  The  predictions  of  the  performance 
curves  reveal  that  a  larger  number  of  flow-channel  leads  to 
better  performance  for  the  effective  distribution  and  utiliza¬ 
tion  of  reacting  gas. 


3.3.  Effect  of  gas -diffusion  layer  porosity 

The  effect  of  gas-diffusion  layer  porosity  in  the  range 
0.3  <  s  <  0.6  on  current  density  distribution  at  the  interface 
between  this  layer  and  the  catalyst  layer  along  v-direction 
of  the  models  with  N  =  2,  and  X  =  0.5  is  shown  in  Fig.  11. 
An  increase  in  porosity  enhances  the  diffusion  transport 
of  the  reactant  gas  through  the  porous  layers  and,  in  turn, 
increases  the  local  values  of  the  current  density  with  a  spec¬ 
ified  overpotential.  This  strong  transport  rate  of  the  reactant 
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Fig.  7.  Comparison  of  /  vs.  V  performance  for  flow  channel  models  with  various  channel  width  fractions. 
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Fig.  8.  Effect  of  surface  overpotential  on  current  density  distributions  along  v-direction  of  models  with  X  =  0.5  and  channel  number:  (a)  TV  =  4;  (b) 
N  =  2;  (c)  N  =  1. 


gas  in  the  across-the-cell  direction  does  result,  however,  in  a 
relatively  non-uniformity  in  the  v-distribution  of  the  current 
density.  Obviously,  cell  performance  depends  on  the  GDL 
porosity.  In  Fig.  12  the  I-V  curves  demonstrate  that  a  better 
performance  can  be  obtained  by  using  a  gas-diffusion  layer 
of  higher  porosity.  The  remarkable  drop  in  cell  potential 
is  caused  by  a  mass  transfer  or  concentration  loss,  which 


is  a  consequence  of  shortage  of  the  reactant  gas  at  high 
current  density  in  which  more  reactant  gas  is  required  for 
fast  reaction.  A  gas-diffusion  layer  of  higher  porosity  has 
an  ability  of  stronger  diffusion  transport,  which  is  beneficial 
in  that  it  supplements  the  reactant  gas  to  the  catalyst  layer 
and  thus  shifts  the  occurrence  of  the  performance  drop  to  a 
higher  value  of  the  current  density  (Fig.  12). 
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Fig.  9.  Correlation  of  local  maximum  of  current  density  and  overpotential  for  channel  models  of  A.  =  0.5  and  N  =  1,  2,  and  4. 
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Fig.  10.  /  vs.  V  performance  curves  for  flow  channel  models  of  X  =  0.5  and  N  =  1,  2,  and  4. 
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Fig.  11.  Effects  of  gas-diffusion  layer  porosity  on  current  density  distributions  along  x-direction  of  models  for  N  =  2,  X  =  0.5  and  various  overpotential: 
(a)  s i  =  0.3;  (b)  s\  =  0.4;  (c)  s\  =  0.5;  (d)  £\  =  0.6. 
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Fig.  11.  ( Continued ). 
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Fig.  12.  Comparison  of  cell  performance  curves  for  models  with  porosity  in  range  0.3  <  e\  <  0.6. 


3.4.  Criterion  for  justification  of  one -dimensional  model 
analysis 

From  the  current  density  results  discussed  above,  it  is  re¬ 
vealed  that  uniform  distribution  of  the  current  density  can 
be  attained  at  a  low  overpotential.  To  facilitate  the  analy¬ 
sis,  a  parameter  is  defined  as  the  relative  deviation  of  local 
current  density  along  the  v-direction  (/max  —  7m in)//max,  to 
judge  of  the  uniformity  of  the  current  density  distribution. 
/max  (=/o)  and  /min  denote  the  maximum  and  minimum  val¬ 
ues.  For  (/max  —  /min) / /max  <  10%,  the  current  density  dis¬ 


tribution  can  be  regarded  as  a  weak  function  of  v.  Therefore, 
the  v-dependence  can  be  neglected  and  an  one-dimensional 
analysis  can  be  regarded  as  appropriate;  otherwise,  at  least 
two-dimensional  analysis  should  be  considered.  At  a  higher 
overpotential,  more  reactant  gas  is  required  for  the  relatively 
strong  and  fast  reaction,  and  therefore,  the  reactant  gas  from 
the  flow  channel  may  not  have  sufficient  time  to  distribute 
uniformly  over  the  width  of  the  cell.  With  non-uniform  dis¬ 
tribution  of  the  reactant  mass  fraction,  the  resultant  current 
density  distribution  also  has  a  noticeable  non-uniformity. 
Based  on  the  computations  with  the  developed  model,  the 
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Fig.  13.  Critical  overpotentials  for  10%  deviation  in  uniformity  of  current  density  distribution. 


values  of  overpotential  for  (7max  —  /min)/7max  =  10%  are 
plotted  in  Fig.  13.  It  is  observed  that  the  critical  overpoten¬ 
tial  is  increased  with  an  increase  in  either  A,  A,  or  s.  A  cor¬ 
relation  of  t)cy  =  Oct  (K  N,  s )  can  be  fitted  with  the  present 
predictions.  Finally,  this  gives: 

r]cr  =  o.i65A00935V0-2033e0-0856  (12) 

The  calculated  values  of  rj cr  by  Eq.  (12)  with  the  parameters 
considered  in  the  present  analysis  are  of  only  less  than  3% 
deviations  from  the  numerical  predictions.  With  a  set  of  the 
given  data  of  A,  N,  and  s,  the  critical  value,  77  Cr,  can  be 
evaluated  from  Eq.  (12).  For  a  value  of  overpotential  ij  <  rj cr, 
the  one-dimensional  analysis  can  be  accepted  for  only  10% 
or  less  non-uniformity  in  prediction  of  the  v-distribution  on 
of  the  current  density.  This  correlation  may  be  employed 


as  a  guide  for  justification  of  the  appropriateness  of  the 
one-dimensional  approximation  in  the  analysis. 

4.  Concluding  remarks 

A  theoretical  model  has  been  developed  to  investigate 
the  effects  of  flow  distributor  and  diffusion  layer  morphol¬ 
ogy  on  the  transport  phenomena  of  reactant  gas  from  flow 
channels  to  gas  diffuser  and  the  effects  on  performance  of  a 
PEMFC.  By  applying  this  two-dimensional  mass  transport 
model,  the  effects  of  overpotential,  channel  width  fraction 
in  the  cross-sectional  view  of  the  anode  flow  distributor, 
number  of  flow  channels,  and  porosity  of  the  gas-diffusion 
layer/membrane  on  the  resultant  current  density  are 
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investigated.  The  cell  performance  is  evaluated  by  the 
calculated  performance  curves.  Based  on  the  results,  the 
following  conclusions  can  be  drawn: 

(i)  Overpotential  can  affect  significantly  both  the  level  and 
the  uniformity  of  the  local  current  density  distribution. 
A  low  overpotential  leads  to  a  more  uniform  current 
density  distribution.  Based  on  the  present  results  at  var¬ 
ious  overpotentials,  a  correlation  is  proposed  to  serve 
as  a  criterion  for  justification  of  the  validity  of  a  sim¬ 
plified  one-dimensional  model. 

(ii)  With  a  fixed  dimension  of  the  flow  distributor  or  bipolar 
plate,  a  larger  channel  width  fraction  results  in  a  better 
cell  performance.  For  a  fixed  channel  width  fraction, 
an  increase  in  the  number  of  flow-channel  results  in  a 
more  effective  gas  distribution,  a  more  uniform  local 
current  density  and,  in  turn,  a  better  performance. 

(iii)  The  morphology  of  the  gas-diffusion  layer  in  the 
present  half-cell  model  is  characterized  by  the  poros¬ 
ity.  With  a  porosity  in  the  range  0.3-0. 6,  a  better  cell 
performance  can  be  obtained  with  higher  porosity  of 
the  porous  layer. 

(iv)  The  present  results  provide  a  useful  guide  for  the  de¬ 
sign  of  a  flow  distributor  and  a  gas-diffusion  layer  via 
a  better  understanding  of  the  effects  of  the  parame¬ 
ters  pertaining  to  the  geometry  and  the  morphology. 


More  profound  investigations  are  needed  to  include  the 
electrochemical  effects  for  a  more  precise  prediction  of 
cell  performance. 
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